Stochastic Resonance in an Overdamped Monostable System. 
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We show, the SR can appear in monostable overdamped systems driven by additive mix of peri- 
odical signal and white Gaussian noise. It can be observed as non-monotonic dependence of SNR on 
the input noise intensity. In this sense it is similar to classical SR observed in overdamped bistable 

■ systems with potential barrier. 

PACS numbers: 05.40.-a, 02.50.Ey 

as: 

■ The dynamics of nonlinear periodically driven stochastic systems has been attracted great attention during the 
last decades . The interest in these systems is much stimulated by phenomenon known as stochastic resonance (SR), 
where noise plays a constructive role p], This phenomenon can be defined as enhancement of sensitivity of a 
nonlinear system to external periodical forcing. Nowadays the SR has been found and studied in a different physical, 

i— i . chemical and biological systems [E 13, IE IE S IE IE EE El • The enhancement of sensitivity is usually understood as 
nonmonotonic dependence of the signal-to-noise ratio (SNR) or signal power amplification (SPA) at the output of 
^ , nonlinear system as a functions of the input noise intensity. Accordingly, the phenomenon of SR displays itself, when 

S SNR and SPA reach maximum at some value of noise intensity and then decrease with further growth of fluctuations. 
The effect of SR was observed in various monostable systems: with multiplicative noise ||, signal array [T^, 
underdamped [l3j], higher harmonics [HI]. The present paper is dedicated to the investigation of monostable systems 
described by the overdamped Langevin equation with an additive noise and an additive driving signal 
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where s(t) — Acos(uJot) is the input driving signal, is the input white Gaussian noise: < >= 0, < £(t)£(i 
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t) >= 2qS(r), 2q is the noise intensity, Q(x) is potential field describing the system and x(t) is the output random 
process. 

The canonical example of SR was observed and studied in overdamped syste m (ffl) for the bistable potential profile 
$(x) with single potential barrier separating the metastable states in Refs. pi 13. [Tof/This result can be generalized for 
multistable potentials with arbitrary number of barriers. The value of additive noise intensity for which SNR reaches 
the maximum was revealed to be about the height of the potential barrier. In other words, the presence of potential 
barrier(s) has been considered as necessary condition for arising of SR in overdamped systems with additive noise fl}. 
j It is well known also that in bistable (multistable) systems the non-monotonic dependence of SNR is accompanied by 
^) . the similar non-monotonic dependence of SPA on noise intensity. 

Recently, in Ref. [ijj it was shown that some special kind of SR can appear in overdamped monostable systems 
qq , ([1]), where there is no any barriers in potential profile $(x). This kind of SR is different from classical one, because 
the non-monotonic dependence on the noise intensity is observed only for SPA. While the SNR was shown to be 
monotonically decreasing function of the noise intensity regardless to various non-monotonic dependencies of SPA 

o ■ M- 

In the present Letter we demonstrate that the non-monotonic dependence of SNR, similar to that for bistable 
systems, can be observed also in monostable overdamped systems. This result implies the presence of SR in monostable 
overdamped system with additive noise and signal. For this case the maximum of SNR for non-zero noise level has 
never been observed before. On the other hand, this situation is sufficiently general to be achieved in a great vaiety 
of physical, chemical, and biological systems. 

The function of SNR is obtained analytically. The analysis of the SNR and SPA shows that this phenomenon of 
SR is new and it has different properties comparing to the classical case appearing in the systems with barriers. 

Consider the following piece-wise linear monostable potential (See Fig. 1) 

1 j I k 2 (\x\-L) + k 1 L, \x\>L. [Z) 
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This potential profile is monostable and has two parameters specifying the slope of potential profile wells: ki 
describes the slope near the minimum at \x\ < L and ki for |x| > L. The both values k\ and ki are always positive 
providing the monostability (single minimum) of potential $(x), while k\ can be grater than ki and vice versa. 

For derivation of the power spectrum density (PSD) of the output signal x(t) the linear response theory (LRT) is 
used assuming that the magnitude A of driving signal s(t) is small enough: A <C minjfc!^}. In accordance with 
LRT (See for example Ref. p]]) the PSD of the output process x(t) reads 

S x (u) = ( w ) + °- (6(lu - luo) + 8{uj + wo)) . (3) 

The function (w) provides the noise platform and the other term is the output signal with amplitude a. Therefore 
the SNR is defined as follows [l[ 

a 2 

R =^F) ■ ( 4 ) 
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The function Sx (w) is the PSD of the unperturbed system (Q]) under s(t) — 0, which is defined as the Fourier 
transform of the appropriate unperturbed autocorrelation function 

oo 

Si Q \u>) = yjKi%]coa{ U T)dT, (5) 



where Kx [t] — (x(t)x(t + r)). In the above expression we have taken into account that JC^Ir] is even function. 
According to the LRT, the amplitude of output signal is 

a = A\x(iuJ )\, (6) 
where x(* w ) is the susceptibility of the system. Therefore SPA of input signal s(t) reads 

V = ^ = \x(^o)\ 2 - (7) 
The susceptibility is the Fourier transform of linear response function h(r) 

oo 

X M= / h^e-^dr, (8) 



while the linear response function can be expressed in terms of correlation function of unperturbed system in accor- 
dance with fluctuation-dissipation theorem 

q dr 
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where 8(t) is the Heaviside function. 

The pro bability density function (PDF) of the unperturbed process W^°'(x, t) satisfies the Fokker-Planck equation 
(FPE) [H] 

?2^ = ^(^<o W)+ (10) 

dx 

with boundary conditions W^°\±oo,t) = 0. Since we consider $(x) — > oo under x — > ±00, the stationary PDF will 
be established in the system with time 

W g\x) = Nexp(-^-), (11) 



q 

where N is the normalization factor. Therefore to find autocorrelation function 



kW[t]= / x W^(x )dx / xW^(x \x,r)dx. (12) 



it is necessary to obtain the transition probability density (xo\x,t), which is the solution of FPE (fT0|) with the 
initial conditions W^°'(x, 0) = S(x — xq). 

For the real physical system the integration in Eq. (|8|) can be not from —00 but from 0, because linear response 
function, according to Eq. ©, exists only if t > 0. Therefore susceptibility x(* w ) m Eq. © can be considered as the 
Laplace transform of the linear response function 

/>oo 

X (p) = Hp) = / h(r)e-P T dr, 
Jo 

where p — iui is the Laplace variable. On the other hand, we can find Laplace transform of the linear response function 
by integrating the expression @ and finally we obtain 



xH^(if(°'[o]-M<°»N), (13) 

where [0] is the correlation function at r = 



K^[0] = (x 2 )+m 2 st , (14) 

which is expressed in terms of variance and mean value of the stationary distribution (fTTj) . The PSD §5§ also can be 
written as the real part of the Laplace transform [p] 

^°)H = iRe{/f(°)H}. (15) 

In the present paper the Laplace transform of autocorrelation function is obtained for monostable potential ([2]) . The 
Laplace transform method for solution of the FPE is described in Refs. 0, [2(J HH [HJ • In particular, in Ref. [20j the 
exact Laplace transform of transition probability density is obtained for piece-wise linear potential profile consisting 
of an arbitrary number of linear parts. Using this approach we obtain the following exact expression for Laplace 
transform of unperturbed autocorrelation function 
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where coefficients Ai are as follows 



xl>] = 7 + ^-. (16) 



q 2 + 2hS + h 2 S 2 - (h 2 + 2h + 2A) S 3 



2pPS 2 1 - XS 



An = 



4qPXfi{l-4S) 



A x = - 



8q0(l-6)((l-\)al+a$) 71 
(1 - Ai) 



A 2 = 



2g (2A/i (2 - 8) - ASfi + SX{2- /x) 71) 
(T^A5) : 



A, = - 



4 g7 i ((2 - 8) a\ + (2 + (2 - fj) 72 )) 
<5 2 (1- A5)(l + 72 ) 



4g7io| (2 - 35) (1 - A) 
8 2 (1 - A<J) (1 + 72) ' 



2g^ ((2 - 5) A + 28 (2 - A) 72 + 25 (1 - A)) 
5 2 (1 - Atf) (1 + 72) ' 



gTj (2 - m) (^A (1 + 72) - 6 + 28 - 2 72 ) 
(38 2 (1 - XS) (1 + 12 ) 



A 7 = 



2qn (1 + 8 (2 - A) - 2<5 2 A) 72 
/^ 2 (l-A<5)(l + 72) : 



As = 



2g/i (3 + 4A - 2<S 2 A - 5 (2 + 3A)) 
/3<5 2 (l-A<5)(l + 7 2 ) 



A, 



/? 2 <5 3 (1 + 72) ' 



B = P 2 (1 + 72 f ((2 - m) 72 - M (1 - $ + Sj2)) , 



here h = hL/q, ( 3 = k 2 /2pq, 71 = i/l + Apg/ k'j, 72 = V 1 + 4 P<?/ fc 2> = fe/fci, A = 1 - e h , /i = 1 - e^ 1 , 
ni = ^1 - fi - \j\J\ — A, «2 = VI — M _ \/i — A- 




FIG. 2: The SNR of the system JIJ with monostable potential (|2| versus dimensionless noise intensity q/H for the amplitude 
of input signal A — 0.3 and frequency luq — 0.1. The Curve 1 corresponds to fci = 25, = 3, L — 0.2 and Curve 2 corresponds 
to fci = 20, fc 2 = 1, and L = 0.1. 



5 



With the above exact expression for autocorrelation function (fT6)) we can obtain the PSD (fT"5|) and the SNR ([4}. In 
Fig. 2 the SNR is plotted versus the dimensionless input noise intensity q/H for different values of L and the slopes 
k%, k 2 . Namely, curve 1 is for L = 0.2, fa = 25 and k 2 = 3, curve 2 is for L = 0.1 and k\ = 20 and k 2 = 1- The 
value H = k\L is the parameter of the potential profile ^ (See the Fig. 1). As one can see from Fig. 2, for small 
and lar ge q , the SNR is decreasing function of q similar to the other cases of monostable potentials considered earlier 
in Ref. [la |. While for the intermediate values of noise intensity 0.1H < q < OAH the non-monotonic behavior of 
the SNR appears. It looks similar to the effect of the SR observed in bistable systems with barrier (See for example 
Refs. [HEl). 

The properties of this SR behavior are varying, depending on the values of parameters of the investigated monostable 
potential profile The effect is observed when ki > 3k 2 . It follows from Fig. 2, that the maximum of SNR is reached 
at noise intensity q* ~ OAH. The difference between the values of SNR in the maximum and in the minimum is 
growing with the difference between fa and fa. In particular, for the Curve 2 in the Fig. 2 we can see the improvement 
in SNR about 10 times with increasing of input noise from q = OAH up to q* = OAH. 

In spite of the similar manifestation, the mechanism of this SR should be different from that in bistable systems, 
where the key parameter of the SR effect is the height of potential barrier. In the considered monostable system there 
is no barrier and the force, which is regular in time, always tends to return the system to the equilibrium point x = 0, 
corresponding to the minimum of potential profile ([2]). The SR appears for the special shape of potential profile, 
which defines the strength of the regular force. Namely, the regular returning force should be much weaker in the area 
located far from equilibrium point comparing to that being near the equilibrium. For the potential ^ the strength 
of the force is changed in the points x — ±L. If this change is large enough fa < 3fa, the SR can be observed as a 
maximum of SNR as a function of the noise intensity. 

Such a shape of potential profile provides the system properties, which are similar to those for excitable systems, 
where the SR also was observed @, [24|. The excitable systems also have only one stable state. Under a small 
perturbation these systems relax quickly to the stable state. While a large (over threshold) perturbation switches 
system to an excited state. The exited state is not stable and it decays to the stable state but after relatively longer 
time. For the potential profile © when k\ ^S> fa the system ((IJ returns to the stable state relatively quickly, if a 
perturbation is less than L. When a perturbation exceeds the threshold value L, the system delays in the region 
> L for a relatively longer time, because the returning force there is weak. On the other hand, the similarities 
between the investigated and excitable systems are only qualitative. The equations describing excitable and threshold 
systems are different from Eq. (JT|). 

n' 
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FIG. 3: The SPA of the system |T} with monostable potential ([2]) versus dimensionless noise intensity q/H for the amplitude 
of input signal A = 0.3 and frequency u>o = 0.1. The Curve 1 corresponds to fa = 25, fa = 3, L = 0.2 and Curve 2 corresponds 
to fa = 20, fa = 1, and L = 0.1. 

The properties of the SR effect observed here for the monostable potential (J2J) has another important difference 
from the SR in bistable systems. In accordance with analysis carried out by various authors, in bistable (and multi- 
stable) systems the non-monotonic behavior of the SNR with the input noise is accompanied by similar behavior of 
the SPA. Using the Laplace transform of the autocorrelation function (fT6|) . we can obtain the SPA for the investigated 
monostable system. The plot of SPA as a function of the input noise intensity is shown in Fig. 3. The parameters of 
the system for the Curves 1 and 2 in the Fig. 3 are the same as those in Fig. 2, where the plots of SNR are presented. 
One can see, the non-monotonic dependence of SNR corresponds to the monotonically growing behavior of SPA as a 
function of noise intensity. 

This diversity in the properties for the SR in bistable and monostable systems confirms the assumption about 
different mechanisms for the SR effect in these systems. Therefore we can conclude that the SR effect presented here 
is new and needs further detailed analysis. 

The authors acknowledge fruitful discussions with Bernardo Spagnolo. This work is supported by Russian Founda- 




tion for Basic Research (project 08-02-01259). 



[1] V.S.Anishchenko, Nonlinear dynamics of chaotic and stochastic systems. - Berlin: Springer 2002. 
[2] L.Gammaitoni, P.Hanggi, P.Jung and F.Marchesoni, Reviews of Modern Physics 70 (1), 223 (1998). 
[3] M.I.Dykman, T.Horita, J.Ross, J.Chem.Phys. 103, 966 (1995). 

[4] P.E.Greenwood, L.M.Ward, D.F.Russell, A.Neiman and F.Moss, Phys.Rev.Lett. 84 (20), 4773 (2000). 
[5] F.Guo, Y.R.Zhou, S.Q.Jiang and T.X.Gu, J.Phys.A: Math.Gen 39, 13861 (2006). 
[6] P.Jung, Chem. Phys. Chem 3, 285 (2002). 

[7] S.Mitaim and B.Kosko, Proceedings of the IEEE 86 (11), 2152 (1998). 
[8] J.M.G.Vilar, R.V.Sole and J.M.Rubf, Phys.Rev.E 59 (5), 5920 (1999). 
[9] S.Zozor, P-O.Amblard, IEEE Transaction on signal processing 51, 3177 (2000). 
[10] D. Valenti, A. Fiasconaro, B. Spagnolo, Physica A 331, 477 (2004). 

[11] B. Spagnolo, M.Cirone, A. La Barbera, F. de Pasquale, J. Phys. Cond. Matt. 14, 2247 (2002). 
[12] J.F.Lindner, B.J.Breen, M.E.Wills, A.R.Bulsara and W.L.Ditto, Phys.Rev.E 63, 051107 (2001). 
[13] N.G. Stocks, N.D. Stein and P.V.E.McClintock, J.Phys.A:Math.Gen. 26, L385 (1993). 
[14] A.N.Grigorenko, S.I.Nikitin and G.V.Roschepkin, Phys.Rev.E 56 (5), R4907 (1997). 
[15] B. McNamara and R. Wiesenfeld, Phys. Rev. E 39, 4854 (1989) 

[16] N.V.Agudov and A.V.Krichigin, International Journal of Bifurcation and Chaos 18 (2008). 
[17] M.Evstigneev, P.Reimann, V.Pankov and R.H.Prince, Europhys.Lett. 65 (1), 7 (2004). 

[18] HRisken, The Fokker-Planck Equation. Methods of Solution and Applications. - Berlin: Springer- Verlag 1989. 

[19] J.D.Atkinson and T.K.Caughey, Int.J.Non-Linear Mech 3, 137 (1968). 

[20] N.V. Agudov and A.N. Malakhov, Radiophysics and Quantum Electronics, 36, 97 (1993). 

[21] A.N. Malakhov and A.L.Pankratov, Physica A 299, 109 (1996). 

[22] V.Privman and H.L.Frisch, J. Chem. Phys 94, 8216 (1991). 

[23] A. Longtin, J. Stat. Phys. 70, 309 (1993). 

[24] A. Longtin and DR. Chialvo, Phys. Rev. Lett. 81, 4012 (1998). 



